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The time evolution of cohective modes in an expanding ultarelativistic and (effectively) 
Abelian plasma is studied in the hard-loop approximation semi-analytically by means of 
integro-differential equations. A previous treatment is generalized to arbitrary orientation of 
' wave vectors with respect to the direction of anisotropy and thus to a fully 3+1 dimensional 

^SJ I situation. Moreover, initial fluctuations are allowed in both gauge fields and currents, which 

^ • is necessary in the case of (stable) longitudinal modes. For unstable (Weibel) modes, this 

^ I generalization of initial conditions reduces drastically the lower bound on the delay in the 

onset of growth that was found previously by considering only collective gauge fields as 
seeds. This makes it appear much more likely that non-Abelian plasma instabilities seeded 
by small initial rapidity fiuctuations could play an important role in the early stage of heavy- 
ion collisions at LHC energies. 
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INTRODUCTION 



(— I ' Fits of hydrodynamical models to the experimental results at the Relativistic Heavy Ion Collider 

(RHIC) are often interpreted as an indication of an extremely fast thermalization, or at least 
! isotropization, of the quark-gluon plasma that is assumed to have formed with initial temperatures 

^> ' significantly above the deconfinement temperature. The inferred thermalization time of < 1 fm/c is 

so short that is hard to understand from a perturbative framework such as the (original) bottom-up 
, thermalization scenario 0-0] • Together with the low inferred value for the specific shear viscosity, 

•/^ I it seems to clearly favor strong-coupling approaches, in particular those based on gauge-gravity 

(N ■ duality f?\. 

However, as pointed out first by Ref. Q], a weak-coupling approach has to take into account 
the inevitable presence of non-Abelian (chromo- Weibel) plasma instabilities [o-JJ], which pro- 



duce nonperturbatively lar ge gaug e fields that may lead to important, qualitative modifications of 



bottom- up thermalization 13Hl7t|. Plasma instabilities and associated turbulent phenomena may 
I even be responsible for a strong reduction of the effective shear viscosity [18|]. A complete scenario, 

' even of just the early stages of the evolution, is still missing, but would clearly be needed to decide, 

firstly, whether the thermalization of the quark-gluon plasma at RHIC is indeed a strong-coupling 
phenomenon, and secondly, what to expect if the early stage of heavy-ion collisions at the higher 
energies at the Large Hadron Collider (LHC) is probing a regime where weak-coupling approaches 
based on (resummed) perturbative quantum chromodynamics (QCD) become relevant. 

At sufficiently weak coupling, the collective dynamics of a non-Abelian plasma can be described 
to leading order by an effective field theory produced by integrating out the hard modes corre- 
sponding to real (as opposed to virtual) plasma constituents. For anisotropic plasmas, the resulting 
hard- loop effective theory [l^ is a generalization of the well-studied hard-thermal- loop effective 
theory Q-H- 

The corresponding effective field equations are nonlocal, but can be made local at 



the expense of introducing a continuous set of auxiliary fields j23l | which arise naturally when solv- 



ing gauge covariant Boltzmann-Vlasov equations 2442 7l|. In the hard loop approximation, these 
auxiliary fields depend on the velocity vector of the hard particles whose hard momentum scale is 
integrated out. 

The instabilities in a stationary homogeneous plasma with momentum-space anisotropy have 
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been studied in this approximation in 0, 2^31] for the case of weak fields, where the dynamics 
is effectively Abelian. In the Abelian case, instabilities grow exponentially until they are large 
enough to modify the distribution of the hard particles, typically giving rise to fast isotropization. 
In the non-Abelian case, nonlinear self-interactions of the collective fields may hinder this evolu- 
tion, which can only be studied by a numerical (real-time lattice) treatment. The first numerical 
simulations of non-Abelian plasma instabilities have only considered modes that are constant in 
the directions transverse to the direction of momentum anisotropy |32l] . In this situation there is 
a short stagnation of the exponential growth when the non-Abelian regime is entered, after which 
localized Abelianization occurs together with continued exponential growth. However, in the more 
generic case of fully 3+1-dimensional evolution [33, 134] such local Abelianization appears to get 
largely destroyed. Instead, a turbulent cascade \2M\ is formed where the growth is reduced to 
a linear one (see however 36| and 37|]). This was also found in more extensive simulations with 
stronger anisotropies in 38|] , although only when initial fields were already nonperturbatively large. 

In Ref. [3^ the hard-loop effective theory for stationary anisotropic plasmas was extended to 
the case of a boost-invariant, longitudinally expanding distribution of plasma particles, which is 
closer to the actual situation in the earliest stage of heavy-ion collisions. In the expanding case, 
more and more modes become unstable as the expansion increases the momentum anisotropy, 
while the growth rate of each unstable mode decreases with the density of the plasma. In [401] 
first numerical results for the non-Abelian evolution were obtained, albeit still restricted to the 
effectively 1-l-l-dimensional situation of modes that are constant in transverse directions. In this 
case, continued approximately exponential growth was observed (albeit only exponential in the 
square root of proper time due to the linearly decreasing density of the plasma). 

In the expanding case, even the time evolution in the weak-field regime is nontrivial and can- 
not be given in closed form. However, in Ref. [s^ integro-differential equations were obtained 
which made it possible to study the time evolution semi-analytically in the case of effectively 
H-l-dimensional dynamics. As a preparation of fully 3+1-dimensional real-time lattice studies 
of non-Abelian plasma instabilities in expanding plasmas, we shall generalize the semi-analytical 
treatment of Ref. to arbitrary wave vectors of the collective modes. 

At the same time, we shall generalize to arbitrary initial conditions, allowing for seeds both in 
the collective gauge fields (the case considered in 0, Q) and also in the auxihary fields corre- 
sponding more directly to colored fluctuations in the hard particle distribution. While in the case 
of stationary anisotropic plasmas little difference was found between the two possibilities for intro- 
ducing seed fields for instabilities, in the expanding case we find that this generalization reduces 



dramatically the uncomfortably long delay of the onset of growth that was observed in [39|, |40| ■ 
As will be discussed in more detail below, this in fact reverts some of the negative conclusions of 
Ref. MM 

concerning the possible role of plasma instabilities in heavy-ion collisions for explain- 
ing very early thermalization. While at RHIC energies with parameters matched from color glass 
condensate scenarios 4]|) there is still somewhat too little room for plasma instabilities that grow 
from small initial rapidity fluctuations, the situation at LHC energies, which appeared somewhat 



marginal in Ref. [39|, |40(], is now much more favorable with regard to an important role of such 



plasma instabilities if the quark-gluon plasma to be produced at the LHC is described by the 
weak-coupling physics underlying the hard-loop-resummed treatment. 



II. HARD-LOOP EFFECTIVE FIELD EQUATIONS WITH ANISOTROPICALLY 

EXPANDING BACKGROUND 

A (sufliciently small) gauge coupling g introduces a hierarchy of smaller momentum scales below 
the scale of "hard" momenta ]p] = of ultrarelativistic plasma constituents. The "soft" scale 
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~ g\/y\p\, where / is the typical hard particle occupation number, is associated with different 
screening phenomena and the various branches of plasmon propagation. To leading order they are 
described by hard (thermal) loop effective theories. When /(p) is anisotropic, the soft scale is also 
the domain of plasma instabilities, which constitute the dominant nonequilibrium effects at weak 
coupling: the associated rates are parametrically larger than any of the scattering processes among 
plasma particles. 

As long as the amplitude of the gauge fields A <C \/7 |p| , the evolution of the plasma instabilities 
is essentially Abelian and can be studied by a perturbative linear response analysis. For a stationary 
anisotropic plasma, the evolution is simply exponential in time until non-Abelian self-interactions 
might hinder further growth when A > y/f |p| and thereby delay the isotropization coming from the 
backreaction of the collective fields on the distribution of hard plasma particles. In an expanding 
plasma, the Abelian (weak field) regime is complicated by the counterplay of increasing anisotropy, 
which favors the appearance of plasma instabilities, and dilution of hard particle densities as well 
as energy densities in soft collective fields. 

In the following we shall first recapitulate the hard-loop effective field equations for an anisotrop- 
ically expanding non-Abelian plasma as introduced in Refs. 39, \^^, and later specialize to the 
effectively Abelian weak-field regime. In view of the numerical simulations that were carried out 
to study the non-Abelian, strong-field regime, the Abelian evolution can be expected to provide 
an upper limit on the strength of plasma instabilities that are seeded by small initial fluctuations. 



A. Boost invariant expanding background 

We assume a color-neutral background distribution of plasma particles /o(p, x, t) which satisfies 

vdfo{p,^,t) = 0, ^;'^=p7/. (2.1) 

This is satisfied trivially in a stationary homogeneous plasma with arbitrary momentum anisotropy. 
In order to describe the earliest stage of heavy-ion collisions in the limit of large nuclei, we consider 
a plasma that expands in one spatial direction (the beam axis). Requiring boost invariance [i^ ] 
and isotropy in the transverse directions, we are led to [H, |3] 

Mp,x) = Mp±,p',z,t) = fo{p±,p",T) (2.2) 

where the transformed longitudinal momentum is 



p'' = 7(/-/5/), f^ = z/t, ^ = t/T, T = Vt^- z\ (2.3) 



with = \J p\ + (p^)^ for ultrarelativistic (massless) particles. 

Introducing the comoving coordinates of proper time r and spacetime rapidity rj through 

t = Tcoshr], (3 = tanhr], 

z = Tsinhr], 7 = coshry, (2-4) 

we are led to a coordinate system with nontrivial metric 

ds^ = gai3dx°'dx^ = dr"^ - d-x]_ - r'^drj'^, (2.5) 

where indices from the beginning of the Greek alphabet will be reserved for the new coordinates. 
Latin indices will from now on only refer to the two transverse spatial directions. Despite the 
nontrivial metric, we shall only use ordinary derivatives when writing for instance Da = da — 
ig[Aa,-] for gauge-covariant derivatives. 
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In addition to space-time rapidity 77, we also introduce momentum space rapidity y for the 
massless particles according to 

= p_L (cosh y, cos i;^, sine/), sinhy). (2-6) 
In comoving coordinates, we then have 



= cosh. T] — sivthr] = p^cos]i{y — rf), (2-7) 
p^ = —p-qlr^ = {p^ cosh T] — p^ sinh rj) /r 

= p'Vt = [p^smh{y - ri)]/T. (2.8) 

Eq. (|2.ip . where the space-time derivatives act at fixed p_L and p^, becomes 

ip" ■ 5a)/o = (2.9) 

with fixed p^ and thus fixed p±,y,(j) as opposed to fixed p". This is solved by any function of the 
form /o(p,x, t) = fo{p±,Pr]), taking into account that pr^ depends on r and ry according to ()2.8p . 
We choose 



fo{p,x) = /iso ( +PVr^so) = /iso ( \IpI + I (2.10) 

'iso 



which corresponds to a locally isotropic distribution on the hypersurface r = Tiso with increasingly 
oblate momentum space anisotropy at r > Tiso but prolate anisotropy for r < Tiso-^ We shall 
also have to consider a lowest value of proper time, tq, where a particle description of the plasma 
constituents begins to make sense. Depending on whether the parameter riso is smaller or larger 
than To we shall consider a plasma that starts with oblate or prolate momentum distribution. 
The particle distribution function (12.10p has the same form as the one used in Refs. 2^, 2^, H, 



341 ] ■ but the anisotropy parameter ^ introduced therein is now space-time dependent according to 

e(r) = (r/Tiso)' - 1, (2.11) 



and the normalization factor N(^) of Ref. [29|, |32| . |34| | is unity. (The anisotropy parameter 6 used 
in Ref. 171] is related to ^ by ^ ~ for large anisotropics.) The behavior ^ ~ r'^ at large r 
corresponds to having a free-streaming background distribution. In a more realistic collisional 
plasma, ^ will have to grow slower than this. In the first stage of the original bottom-up scenario 



[13| . ignoring plasma instabilities, one would have had ~ r^/^. In Ref. 14] it was argued that 



plasma instabilities reduce the exponent to ^ ~ t^^^, whereas Ref. 17] recently presented arguments 
in favor of ^ ~ r^/^. All these scenarios typically consider ^ ^> 1, so below we shall mostly 
concentrate on the case Tiso < tq and thus high anisotropy for all r > tq, but in the simplified case 
of a collisionless free-streaming expansion. However, we shall also discuss instabilities in prolate 
phases for which we need to set Tiso > ''"o- 



^ Notice that riso is just a parameter of the background distribution and does not refer to the time where isotropiza- 
tion of the plasma eventually occurs through interactions. 
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B. Hard-expanding-loop effective field equations 



In an approximately collisionless plasma, the gauge covariant Boltzmann-Vlasov equations for 
color charge carrying perturbations 6fa have the form 

V ■D6fa{p,^,t) = gv^FrdiP^fo{p,^,t). (2.12) 

In comoving coordinates we write 

V ■ D6f\, = gV''F^^d^^^^fo{p^,p,), (2.13) 

where in place of the light-like vector = /p^ containing a unit 3- vector we introduced the new 
quantity 

= — = cosh(y — ry), cosi^, sinc^, — sinh(y — ry) , (2-14) 
P± \ J 

which is normalized so that it has a unit 2-vector in the transverse plane. 

Eq. (j2.13p can be solved in terms of an auxiliary field Wjs^x; (p, y) which satisfies 



V-DWp\^^y = V''Fp^, (2.15) 

and 



6f{x]p) = -gWp{x;(l),y)d^^^fo{p_L,p^). (2.16) 
Expressed in terms of the auxiliary field W, the induced current in comoving coordinates reads 

cfip 



/dp 



where is a suitably normalized group factor. Now for each {(f),y) (i.e., fixed v) the scale p± 
(related to energy by p^ = p± coshy) can be integrated out. 

With the particular background distribution function (j2.10p we obtain 



where 



2 _ 2. r dpp\, 

(2vr)2 



ml = -gHRl j^fLip) (2.19) 



equals the Debye mass squared at the time Tiso (which we shall often choose smaller than tq so 
that niD is not physically realized by just a convenient mass parameter). 

These equations are closed by the non-Abelian Maxwell equations which in comoving coordinates 
read 

-D^irF"^) = -D^ \rg^nT)g^\T)Fj = f, (2.20) 
r r L 'J 
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where -F^^ = daA^ — dpAa — ig[Aa, Aj^]. To solve them, we adopt the comoving temporal gauge 
A"^ = and introduce canonical conjugate field momenta for the remaining gauge fields according 

to2 

= TdrAi = -TdrA' = -Ui , (2.21) 

and 

n'' = ^drA^ . (2.22) 

In terms of fields and conjugate momenta, the Yang- Mills field equations then read 

rdrU^ = jr, - DiF\ , (2.23) 
T-^drlii = f - DjF^' - Dr,F''' . (2.24) 

In a comoving frame, the longitudinal (chromo-)electric and magnetic fields are given by 

E'^ = m, Br, = Fi2, (2.25) 

but transverse components involve a factor of r, 

= r~^rf, Bi = T'^Frjjeji. (2.26) 

In terms of these, the contribution to the energy density is simply 

£ = £t + £l = £bt+ £et + £b, + £Er. = tr [{Bif + {E'f + (5^)^ + {E^'f] . (2.27) 

However, due to the expansion, the total energy density £ is not conserved, even when the induced 
current ()2.17p is identically zero, 

-^^\j=o = --ST\j=o ■ (2.28) 

In the presence of a plasma of hard particles and thus nonvanishing induced current j we define 
the net energy gain rate by 

d£ 2 

i^Gain = -r + -£T, (2.29) 

which gives the rate of energy transfer from the free-streaming hard particles into the collective 
chromo-fields and which is positive when plasma instabilities are at work. For stable modes -Rcain 
oscillates about zero. 



III. TIME EVOLUTION OF GAUGE FIELDS IN THE WEAK-FIELD REGIME 

In the regime where self interactions of gauge fields cannot be neglected, the Yang-Mills field 
equations and the equations of motion for the VFa(x; ^, y) fields are nonlinear and require numerical, 
real-time lattice evaluation. In the limit of small field amplitudes, these equations are linear and 
can be reduced to ordinary integro-differential equations in proper time for individual modes. 



^ In the following field equations we keep the index position of conjugate field momenta opposite to that of the 
associated fields. 
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A. Solving the W field equations 

In the weak-field regime, the W field equations reduce to 

{V^dr + V'di + V^d^)Wa,{T,x\r^-4>,y) = V^F^p = V^{daAp - d^Ac,). (3.1) 

Because our background distribution function is isotropic in the transverse plane, we can without 
loss of generality restrict to modes which are independent of x*^^, keeping only x*^^ = x (recall 
that the symbol y is already used for momentum rapidity). In temporal axial gauge we thus have 

{V^dr + V'^dr, + cos4>d^)Wo,{T,x,iTA,y) = -V^ drA^ + [d^^A^ - O^A^) 

+ cos(t) {daAi — dxAa) + sm.(f) daA2. (3.2) 

These first-order partial differential equations can be solved by the method of characteristics as 
follows. We introduce a parameter s such that the left-hand side of Eq. (j3.2p is replaced by 



dWo, dWadr dWadrj dW^dx 



ds dr ds drj ds dx ds 

with 



dr 
ds 

? = V^ = ^smHy-rj{s)), (3.5) 
ds T(S) 

dx 
ds 



= cosh{y -r]{s)), (3.4) 



V^ = cos(t>. (3.6) 



Since dr/ds > we can use r in place of s for the purpose of integrating Eq. ()3.2p . Writing 
ds = dr' /V^ {rj{T')) we obtain 

H^«(r,x,r?;,/.,y)-T^,(ro,xo,r?o;</.,y)= / ^t' ^^ ^'fj-' (3.7) 



J TO cosh(y - r/(r')) 
with xq = x{t' = To) and t/q = ^^(t' = tq). The functions ^/(t') and x[t') are solutions of 

^ = itanh(y -r?(r')) (3.8) 

dx{T') ^ COS0 

dr' cosh(i/ — r]{T')) 
with initial conditions 7]{t' = t) = t] and x[t' = t) = x. Eq. (|3.8|) is solved by [sol ] 

r'sinh(y — ij^t')) = Tsinh(y — r]) (3.10) 

or, more explicitly, 

r]' = r]{T') = y — asinh ^— sinh(?/ — rj)^ . (3-11) 

With this solution, Eq. (j3.9p can be integrated, yielding 

x' = x(t') = x + [t' cosh(y — rj') — t cosh(y — r/)] cos (/). (3.12) 
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The W fields, from which the induced current is obtained upon integration over (j) and y ac- 
cording to Eq. (I2.18p . are now given exphcitly by the fohowing "memory integrals" 



Wi - 
W2 - 



dr' 



TO 



sin( 



dr' 



TO 



COS ( 



cosh(?7' — y) 



cosh(77' — y) 



dr' 



dr'Ar, + 



V'd^^A' + d^>A^ 
cosh(?7' — y) 



(3.13) 



where inside the integrals A^ = Aq(t', x', 77'). Note that x' and ij' as well as xq and rjQ appearing 
as arguments in = Wa{TQ, xq, 4>, y) all are functions of the spacetime variables r, r] and x. 



B. Fourier components 

Because of the linearity of the Maxwell (linearized Yang-Mills) equations and the W equations 
in the weak-field (Abelian) regime and their translational invariance in r] and y, we can study the 
time evolution of individual modes obtained by a Fourier decomposition 

dk .^^ r du 



A^{t,x,ii) = j -e'^y —e^-VA^{T-Kv). (3.14) 

With similarly Fourier transformed currents, the equations of motion for A^ and A^ read 

[r'^drTdr + T^'^i''^)A^{T;k,u) = j^{T;k,u) - kuT~'^An{T;k,iy) 

{TdrT~^dr + k'^)Ar^{T-k,u) = jnir; k, 1^) - ki^A^ (t] k, 1^) (3.15) 

with j^(r; A:,i/) and j^{T;k,u) both depending on A^{T';k,u) and yl^(r'; A;, z^), r' < r, but not on 
A^{t' ■,k,v), as we shall see shortly. The equation of motion for the transverse mode A^{T]k,v) is 
decoupled from A^ and and reads 

{T'^drTdr + k"^ + T'"^ v"^) A^ {t ] k, v) = f{T] k, u) (3.16) 

with j^(r; k, v) depending on the history of A^{t'] k, u). 

In order to express the current as functional of the gauge fields, we first note that in Eqs. ()3.13p . 
the partial derivatives dx' and dff' are simply replaced by factors ik and ii^, respectively, while 
the proper-time integrals over the partial time derivative of the gauge fields ^^(r', x', r/') can be 
integrated by parts, yielding 

dT'dr^AUT',x',r,') = f dr' I ^e^'^-'C-')/ ^^'^V {t') q^, X^^^.j^^y) (3.17) 



TO JTO 

dk f du 



ikx^iun^^^^- k, v) - e^'^^o(^)e^'^^«Wi„(ro; k, v) 

Aair'; k,u)\. 



TO 



^ COS (p tanh(y — rj') 

cosh(y — rj') t' 



Inserting Eqs. (|3.13|) with (|3.17|) into the expression for the current, Eq. (j2.18|) and introducing 
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the abbreviations^ 



y = y -V, 

f 

f]' = fi{T') = r]{T') - r] = y - asinh (— sinhy), % = v{to), 

T 



x' = Xir') = (x(r') - x) / cos (j) = \J r'2 + sinh^ y -t cosh y, Xo = x(to), (3.18) 

we find the following results for the Fourier components j"(t; k, v) after performing the integration 
over the (momentum space) angle (j) in terms of Bessel functions Jn, 



m 



D 



irsinhw ^ ,, xt . x Irx -ri/, x xi / x 
^Ji(fcXo)^r,(ro) - -[Jo - J2](A:Xo)^^(ro) 



+ / dr' 



TO , / 1 + li^i^ 



^ ii/rsinhw- ,,\ -i , ,, 

-[3Ji - h]{k^) - [-^0 - J2](/cx')J^'(^') 



Tsinhw/fA;.^ z^rsinhy ^ ,,\ ~ , 
^(^[Jo - J2]{k)(') -^Ji{kt))A^{T') 



(3.19) 



and 



m|) /" dy sinh y 



2r 



r sinh y 



1 + 



sinh^ « \ 2 



i^(r) +e 



r sinh y 



jQ{kxo)K{To) - i'Ji{kxo)^^{TQ) 



+ / dr'- 

Ito Jl-\- '^'^h^ g 



z^r sinh y 



^[J2 - Jom') + ^-^^Jl{kx'))A\T') 



+ - 



r sinh y / iz^r sinh y 



r'2 



Jo(A:x')-A^Ji(A:x))^^(t') 



+ J?(r). 



(3.20) 



As mentioned above, these components of the current only depend on (the history of) and Arj. 
If either k = oi v = 0, the 1 and r] components decouple from each other. (For A; = only the 
terms involving Jo(0) = 1 survive, while for = all terms involving odd powers of sinhy integrate 
to zero.) On the other hand, j'^ is found to be a functional of only A^, 



^2 



m 



D 



dy 



1+ '""r 



(3.21) 



+ / dr'- . 



1[Ji + Jsm') - '-^^^^^[Jo + J2m'))A\r') +ji{r) 



Eqs. (|3.19|) - (|3.2ip generalize the expressions given in 139] for the effectively 1+1-dimensional 
case k = 0, where the wave vector of the collective modes points in the direction of momentum- 
space anisotropy. (As they should, the functional dependences in and become the 
same for k = 0.) 



^ The notation is chosen such that a bar indicates a dependence on y and a prime a dependence on t' {y and r' are 
the two remaining integration variables in Eqs. H3.19p - (|3.22p '). 
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The T component of the current, which is needed only for a check of the Gauss law constraint, 
is again independent of A"^ and given by the following functional of and A^, 



m 



D 



+ 



dy cosh y 

/-^ _l_ sinh^ ^ \ 2 



T sinh y 



Jo{kxo)Ajj{To) - iJi{kxo)A^iTo) 



dr' 



TO 



1 I sinh^ y 
J- H r72 



^^FT Ti/, /N z^Tsinhy _^ , ,,\ ~i , 



rsinh-y /ii^r sinhw ^ , /, /s\ t / 

+ ^( ^2 ^ /o(fcxO-fe^i(fcxO)-4^(TO 



(3.22) 



In all these expressions, ("T"! ^) corresponds to nontrivial initial data for Wa, if any. Up to a 
factor e^^^e^'^v^ the Fourier component jQ(r;fc,i^) is obtained by evaluating the expression for the 
induced current, Eq. (|2.18|) . with VF's of the form 



Wa(To, xo(t), r/o(T); 0, y) = e^'^^^^e^'^^" z.; y - r?o(T)), 



(3.23) 



where the functions xq and r/o are given by Eqs. (I3.12p and (13. lip , respectively, with r' = tq. 
Nontrivial initial values are required whenever the right-hand side of the Gauss law constraint 
at r = To, 



firo) 



To 



zvn^(ro) + m^(To; 



(3.24) 



is nonvanishing. This is naturally the case for modes where the polarization of the gauge fields and 
their momenta (in temporal gauge) is longitudinal with respect to the spatial wave vector. 

Note that already for a single mode there is considerable freedom in the choice of initial con- 
ditions which is parametrized by the functions W^{k,i';(l),y — 7]q) and which will be explored in 
detail below. 



C. Stable and unstable modes 

Before studying the time evolution of the individual Fourier modes through numerical evaluation 
of the integro-differential equations provided by the above expressions, it is useful to recall the case 
of a stationary anisotropic plasma, where a rather complete analysis of stable and unstable modes 
in the regime of weak fields has been carried out in Refs. 0, [2^ . [2^ . This is still relevant for the 
expanding case, which is however complicated by the time dependence of the density of the plasma 
and its anisotropy parameter so that some unstable modes may shut off as the plasma evolves 
towards a higher degree of oblateness, while new ones come into being. Since the growth rate of 
unstable modes also depends on the orientation of the wave vector, a further complication comes 
from the fact that this orientation is time dependent unless the wave vector is strictly parallel 
or orthogonal to the anisotropy direction. When both k ^ and 7^ 0, in comoving Cartesian 
coordinates the wave vector rotates into the transverse plane according to 

k = A;ei-F-e3 (3.25) 
r 

with the angle between k and 63 given by 

^? = arctan— . (3.26) 
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TABLE I: Classification of modes according to Ref. [28|, the corresponding polarization of the electric and 
magnetic field, and the range of instabilities in terms of the angle d between the wave vector k and the 
direction of anisotropy. 



Mode 


E field 


B field 


instabilities, prolate case 


instabilities, oblate case 


a 


k_L E II e2 


k _L B _L 62 


stable 


0<^l< f 


+ 


E±e2 


k_L B II 62 


stable 


stable 




E_Le2 


k±B II 62 


iL < ^ < ZL 
4 ^ - 2 


< I? < f 



When the anisotropy is characterized by only one spatial direction, as is the case for the dis- 
tribution (j2.10p . there are in general three different branches^ of modes. Following Ref. [28], the 
modes with a polarization of the electric field transverse to both wave vector and the direction of 
anisotropy are denoted by the label a. In our case these correspond to the Fourier components 
A^{t; k, u) and j^(t; v), described by Eqs. (|3.16p and (j3.2ip . Such modes are stable in the case of 
a prolate momentum anisotropy, and unstable for all orientations of the wave vector in the case of 
oblate anisotropy, but the growth rate of the latter approaches zero as — )• ^ . In the terminology 
of Ref. ^] , these instabilities are magnetic ones and have been first described by Weibel [i^ . The 
instabilities studied so far for expanding plasmas 39,(43] correspond to such modes in the special 
case 1? = 0. 

When the polarization of the electric field (and of the gauge field in temporal gauge) lies in the 
plane spanned by the wave vector and the direction of anisotropy, there are two modes, labelled 
"+" and "— " in Ref. 28|], according to which one has the larger (+) or smaller (— ) zero-frequency 
mass squared in the stationary anisotropic case. In the isotropic case, the mode corresponds to 
electric Debye screening at zero frequency and longitudinal plasmons above the plasma frequency, 
whereas the "— " mode coincides with the a mode (which in this case has zero screening mass 
to leading order). In the anisotropic case, this degeneracy is lifted as soon as 7^ 0. Now the 
"— " mode is unstable for both oblate and prolate anisotropics, but in each case for only a limited 
range of angles t?, see Table HI 

For generic orientation of the wave vector, the "— " mode involves also longitudinal electric 
fields and thus can correspond to an electric instability according to the classification of Ref. @|. 
However, when either A; = (t? = 0) or 1/ = (t? = ^), the equations for and A^, (j3.15p . 
(j3.19p and ()3.20p . decouple. The mode which is longitudinal with respect to k is then electrically 
screened, whereas the other one is transverse and corresponds to either magnetostatic screening 
(when = ^ in the oblate case, or = in the prolate case), or a magnetic instability (when 
= in the oblate case, or i? = ^ in the prolate case)^. 

In the following we shall study the time evolution of some representative cases of stable and 
unstable modes in the expanding case, for initially oblate as well as prolate anisotropics, with a 
particular view on the dependence on initial conditions. 



* For a given wave vector, one branch may have more than one mode, e.g., a propagating wave and a growing 
unstable mode. 

^ For large anisotropies, the magnetostatic screening mass may even become larger than the electrostatic (Debye) 
mass, e.g. with prolate anisotropy ^ ^ —0.88, as is the case in Fig. 2b of Ref. [28| . whereupon electrostatic Debye 
screening at i9 = becomes part of the "— " branch and thus continuously connected with the magnetic instability 
at = f . 
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IV. NUMERICAL EVALUATION 

The integro-differential equations for the time evolution of individual Fourier modes that we 
have obtained in Sect. IIIIBI for the linear response regime can be solved straightforwardly by 
discretizing the proper time variables r and employing a leap frog algorithm for gauge fields and 
conjugate momenta, Eqs. ()2.2ip - (j2.23p . The memory integrals in the expressions for the induced 
currents also involve integrations over the momentum rapidity variable y that have to be performed 
for each time step between tq and r. 

Since we are interested in the earliest stage of the evolution of plasma instabilities from small 
initial fluctuations, we choose our dimensionful quantities such that tq ~ Qs, with Qg the so-called 



saturation scale of the color-glass-condensate (CGC) framework 4l|, |46|], and a normalization of 



the hard particle distribution function ()2.10p such that at r = tq the hard-gluon density of CGC 



estimates is matched. This involves the so-called gluon liberation factor c 471], which we choose 



as c = 21n2 ~ 1.386 as obtained in approximate analytical calculations by Kovchegov in Ref. 48l |. 
While this value is significantly larger than the first numerical estimates 49, [sO], it is fairly close 
to the most recent numerical result c ~ 1.1 by Lappi 511]. Our choice of c corresponds to a value of 
the squared mass parameter i^'^ in the expression for the induced current of ~ l*285(TQTiso) ^ • 
(See Appendix C of Ref. [io^ for details.) Unless stated otherwise, this value will be used in the 
following numerical calculations. The only remaining free parameter is then Tiso, parametrizing the 
amount of anisotropy at the initial time tq, and we shall consider both initially prolate and oblate 
distributions. 

For later reference we note that if we assume T(^^ ~ ~ 1 and 3 GeV for RHIC and LHC 
experiments, respectively, 1 fm/c corresponds to ~ 5ro for RHIC and ~ 15ro for LHC. 



A. Wave vector parallel to anisotropy direction 

Fourier modes with A; = and u ^ have a wave vector parallel to the spatial direction of 
momentum anisotropy and thus are constant in the plane transverse to the axis of expansion. Such 
modes are stable for prolate anisotropy, while with oblate anisotropy there are magnetic (Weibel) 
instabilities below a certain ((^- or r-dependent) value of v. The case fc = 0, which has been studied 
semi-analytically before in Ref. [39], is particularly interesting since it covers the most unstable 
mode of a plasma with oblate anisotropy. 

Before studying the unstable modes in more detail (eventually also with A; 7^ 0), we begin with 
the stable longitudinal modes, which are complicated by the need for nonvanishing initial induced 



currents. The analysis of Ref. 39|] of the unstable modes at A; = will then be generalized by 



allowing also for nonvanishing initial currents. 



1. Stable (longitudinal plasmon) modes 

Longitudinal modes with k = and nonvanishing A^, H*^ are purely electrical and correspond 
to charge density waves (longitudinal plasmons). Initial conditions only in A^, with zero initial 
H'' and W fields would only yield a trivial, constant solution. For nontrivial solutions, we need 
nonvanishing initial W fields and nonzero initial currents. In order to have nonzero j^{TQ',v) we 
need initial values of the W fields that are odd in y = y — rj. A possible choice is 



W {k = 0, 1/; (/), y - %) = Ci tanh(y - r/o). 



(4.1) 
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FIG. 1: Longitudinal modes with ly = and 1 for Tiso/To = 10. The dotted and dash-dotted lines represent 
the analytic results for the late time behavior, Eq. (jAlSI) . which are in good agreement with the numerical 
data (thick full lines) at sufficiently large times. 



where Ci is a constant. We recall that tjq is given by (|3.1ip as 

Vo = v{^' = 'To) = y — asinh ( — sinh(t/ — ??) ) (4-2) 

so that r]Q = r] at T = tq. The Gauss law constraint p.24p relates the constant Ci to the initial 
value of the longitudinal electric field according to 

27rCim^To f dysinh^y iuIT^ {tq; ly) 



3 {tq;v) = —2 / , . „_,2 = • (4-3) 

The nonzero field (j4.ip thus gives rise to the following contributions jQij) in the integral 
equations (IHT^ and (I22D, 

1 I sinli^ y\'^ Li I sinh^ « 

v 

1 H =2 — - ) A/lH 



where^ 

Recall that f]o = tjq — r], which vanishes at r = tq, so that Jq (''"o! i^) = 



, X /" sinh^ y , . , 

iV(Tiso/T0) = / , ,.J,2 - (4-6) 

' ' Sinn y ^ 



^ See Eq. (|A13|I for an analytic expression for A'^. 
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With A; = 0, the complete expression for the induced current (j3.20p reads 



ISO 



m j-) f dy sinh^ y 



ISO 



dr' i,(.» Uj;(.;.). (4.7) 

TO t'^a/i + '^'"^ y ) 



which has to be solved together with 

dr-drA,,{r;v) = -Tf{T;v). (4.8) 
r 

In Appendix lA 21 the late time behavior of the solutions is derived in terms of Bessel functions 
J2 and 12, see Eq. (jAlSp . In Fig. [1] this is compared with results of a full numerical solution of the 
above integro-differential equation for Arj using t^so/tq = 10 so that the time evolution starts in 
the prolate phase, and v = Q and 1. The late-time (large oblate anisotropy) behavior is reproduced 
very well, with noticeable deviations at earlier times. 



2. Unstable (transverse) modes 

With /c = 0, there are two degenerate transverse modes described by the then coinciding equa- 
tions (j3.19p and p.2ip for the induced currents and j^, which reduce to 

ir' '--f-f^y i-(r-;.)|+j-(r;.). (4.9) 

^2 A _L t2 sinh" y 



ISO 



7>T 



This has to be solved together with 



ili{T-u) = TdrA\ \fi^{T;u) = ~nT;u)-'^A\T;v). (4.10) 



In the notation of Table [U the solutions of these equations correspond to the modes "— " and 
;", which become the same at "i? = 0. This case, which contains magnetic Weibel instabilities 



for oblate anisotropies, was studied already in Ref. 39|], but with vanishing initial W fields and, 
correspondingly, vanishing initial currents, which now is perfectly consistent with the Gauss law 
constraint j'^ oc 11'' = 0. 

In the following numerical evaluations we shall be more general and consider separately the 
initial conditions of (i) a seed electric field with only 11* (tq) 7^ 0, which is the case studied semi- 
analytically in Ref. [39(], (ii) a seed magnetic field with only A^{tq) 7^ 0, which was covered before 
in the real-time lattice calculations of Ref. [13], and (iii) only j*(to) 7^ 0. (In the present linear 
response analysis, the most general case is given by a linear superposition of these possibilities.) 

A nonvanishing initial current io('^o) is provided by any nonzero function Wf{k = 0, z^; y — rjo) 
that is even in its last two arguments. In the following we simply take a constant 

W,°(A: = 0,z.;(/>,y-7?o) = C2, (4.11) 
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FIG. 2: Magnetic (Weibel) instabilities with fc = for an oblate momentum distributions from the beginning 
with Tiso/To = 0.1. In the upper panel transverse magnetic fields are compared for different wave numbers v 
but equal initial conditions n'(To) — 1 and v4*(to) = j'(T"o) = 0. The expected analytical late time behavior, 
Eq. (|Alip . is indicated by the thin light lines. In the lower plot the influence of different initial values is 
studied for modes with v = 30. 

which we found to be also representative of some more complicated possibilities that we have 
studied. Proceeding as above, this determines the function Jo(t) in Eqs. (j3.19p and (j3.2ip with C2 
proportional to the initial values f{To). 

In the upper part of Fig. [2] we consider the case (i) of an initial seed electric field for Tiso/tq = 0.1 
and compare with the analytic result for the late-time asymptotics (jAlip for various values of v. 
This is the analog to Fig. 1 of Ref. [sol but with our larger mass parameter^ and less extreme initial 
anisotropy. Like Ref. [39] we observe a substantial delay in the onset of plasma instabilities which 
is preceded by a decay of collective fields until ~ IOtq, suggesting an uncomfortable suppression of 
Weibel instabilities by the initial strong expansion of the plasma. 

In the lower part of Fig. [21 the dependence of this behavior on initial conditions is displayed for 



As mentioned above, the recent CGC results [5l[ now favor this larger mass parameter, which was also considered 
in Ref. (39I ] but without corresponding plots. 



: \b{t,V)\ 
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FIG. 4: The total energy density Etot for the mode with v — IQ and /'(tq; v) ^ Q of Fig.[3]and its contributions 
from electric {£e) and magnetic fields {£b)- Additionally the gain rate (times tq) is shown. 



mode v = 30. Case (ii) corresponds to using seed magnetic fields instead of seed electric fields, but 
this only increases the delay of plasma instabilities, which is in line with the results of (4o| where 
mixed initial conditions for the fields were considered. Surprisingly enough, with case (iii) which 
corresponds to initial fluctuations in the currents, we find that this delay is very strongly reduced. 

In Fig. [3] the comparison between cases (i) and (iii) is repeated for various values of v, showing 
now the total energy in the collective fields, which confirms the finding of a drastic acceleration of 
the onset of plasma instabilities when there are initial current fluctuations. In Fig. |4] the energy 
density of one of those quickly growing modes {v = 10) is decomposed into magnetic and electric 
contributions together with the gain rate defined in Eq. (12.29p . 

From Fig. [3] one can also easily see the effect of having initial fluctuations in both induced 
currents and (chromo-)electromagnetic flelds. Because everything is linear, the resulting solutions 
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are just linear combinations, and because the onset of exponential growth is so much quicker for 
the part corresponding to initial current fluctuations, such superpositions are dominated over- 
whelmingly by the latter as long as the two components have comparable initial energy densities. 
Since in the physical context of heavy-ion collisions we expect to find fluctuations in all quantities, 
induced currents {W fields) as well as gauge fields, we thus consider case (iii) as being actually 
representative of generic situations. 

The strong reduction of the delay of the onset of Weibel instabilities makes it appear much more 
likely that plasma instabilities could play an important role in the very early dynamics of a quark- 
gluon plasma, at least with LHC energies, where Qs ~ 3 GeV. Choosing tq ~ Qj^ ~ j^fm/c and 
judging from the time it takes that the initial depletion of energy in the fastest mode is reversed 
in Fig. [21 one may set the scale where plasma instabilities kick in to ~ 0.5 fm/c, whereas the 
less generic initial conditions with only seed fields and no currents considered previously in Refs. 
"sil . 40 1 would have given ~ 3 fm/c. (RHIC energies would give values 2-3 times higher.) 



B. Wave vector perpendicular to anisotropy direction 



Another comparatively simple case is provided by a wave vector which is strictly perpendicular 
to the anisotropy direction, i.e. = and k ^ 0. In this case the integro-differential equations 
for A^, and again decouple, but the equations for and AP' are no longer identical, so 
that we have 3 different modes. Now there are two stable modes: the purely electrical mode with 
polarization along the wave vector, described by ^ and the ("a") mode A? ^ where the electric field 
is transverse to both the wave vector and the anisotropy direction and the magnetic field pointing 
in the anisotropy direction. The third mode, j4^, where the electric field points in the anisotropy 
direction while the magnetic field is transverse to both the wave vector and the anisotropy direction 
is stable for oblate anisotropics, but contains a magnetic Weibel instability for the prolate case. 
We shall therefore now consider Tiso > so that we have an initial period of prolate anisotropy 
before the expansion changes that into an oblate one, where none of the 1^ = modes is unstable. 



1. Stable modes 

The purely electrical mode A^ with polarization along the wave vector again requires initial W 
fields in order to satisfy the Gauss law constraint. One of the simplest choices is 

W?ik,i^ = 0;^,y-7jo) = C3Cos(^ (4.12) 

with a constant C3 that is proportional to the initial value ni(ro; k) appearing in the initial Gauss 
law constraint 

ro~f{ro;k) = -ikUi{To;k). (4.13) 

Mode A^, which is transverse to both wave vector and anisotropy direction, does not need 
initial W fields to satisfy the Gauss law constraint. Because we are more interested in the influence 
of different initial conditions on the evolution of plasma instabilities, we shall only consider the 
simplest case of IP{to) / 0. 

In Fig. Owe compare the numerical results for the electric fields corresponding to the two stable 
modes A^ and A^ for k = 0.1rJ~^ and Tiso/tq = 100, with a time range so that both prolate 
and oblate anisotropics are appearing. As expected, only stable oscillatory behavior is found. We 
observe that the purely electrical (plasmon) mode has a smaller frequency than the transverse stable 



mode, which is qualitatively similar to the familiar behavior in the isotropic case 53]. Perhaps 
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FIG. 5: Stable modes for wave vectors perpendicular to the anisotropy direction with k = 0.1 Tq'"'^, initial 
condition Ili{TQ; k) = 1, and Tiso/To = 100. The full line corresponds to the longitudinal plasmon mode , 
the dashed lines to electric and magnetic fields of transverse plasmons. 

more surprising is the rather strong attenuation during the first oscillations which is even stronger 
in the longitudinal plasmon mode. 

2. Weibel instability during prolate phase 

In the time evolution of Arj{T]k,u = 0) we expect to find magnetic Weibel instabilities for 
prolate anisotropies, thus only as long as r < riso- In contrast to Weibel instabilities for oblate 
anisotropics there is now only one unstable mode, namely with electric field transverse to the wave 
vector and pointing along the direction of anisotropy, and the magnetic field transverse to both. 

We again consider initial conditions with and without initial currents. Nonzero initial induced 
currents can be set up most simply by choosing a constant W^{k, u = Q;(f),y — rjo) = C4, which 
we compare with vanishing and nonzero initial electric field, II^{tq) 7^ 0. 

As before we choose Tiso ^ tq to have an extended phase where the free-streaming plasma 
has prolate anisotropy, and we again stop the numerics only after the momentum distribution has 
become oblate. 

Unfortunately the numerical solutions show rather little activity with the mass parameter mo 
extracted from CGC calculations. We therefore choose the much higher value = 1000/(roTiso) 
at first and consider the situation for our standard value only thereafter. Numerical results for 
the higher mass parameter and rjso = 10 tq are displayed in Fig. [6^, which exhibit a pronounced 
instability that shuts off when the degree of prolate anisotropy becomes too small for a given value 
of k. After that point in time the mode decays for about as long as it was growing initially, ending 
in stable oscillations. Dashed lines correspond to nonzero initial values of 11'' (tq) and full lines to 
nonzero initial currents 7^ 0, normalized such that the two different initial conditions have equal 
amplitude in the final oscillations. Since the solutions with nonzero initial currents reach larger 
maximal values we again find, although to a lesser degree, that such initial conditions are more 
efficient seeds for unstable modes. 

In Fig. [6b the energy content in magnetic and electric fields is shown for one of the unstable 
modes, demonstrating that the energy is predominantly in magnetic fields, as expected for a Weibel 
instability. Near the point where the instability stops the electric field changes sign. 
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FIG. 6: Magnetic Weibel instabilities for initially prolate momentum distributions with Tiso/To = 10 for 
different wave numbers k (in units of Tq^) and two different initial conditions (upper panel). The dashed 
lines correspond to Il^{To;k) ^ and the full lines to j'^{To;k) ^0. In the lower panel the total energy 
density and its contributions from electric and magnetic fields are shown for the mode with fc = 10 Tq"^ and 
j''(ro; k) 0. The mass parameter in both plots has been increased to mjy = 1000/(riso''o). 



In Fig. [7] we finally consider our much smaller standard choice mj^TQTiso = 1.285 and again the 
two initial conditions of nonzero initial electric field (upper plot, where W{to) = 1) and nonzero 
initial current (lower plot), normalized so that the final oscillations have equal amplitude to first 
case. Also shown is the gain rate, defined in Eq. (j2.29p . times tq. Notice that the plot is now linear 
instead of logarithmic. In order to observe some instability, we need to consider much stronger 
initial anisotropy and k <^ Tq^. Fig. [7] shows the various energy components for Tjgo = lOOro and 
k = 0.1t-\ 

Since the instability is now a magnetic one, its presence is best judged from the magnetic energy 
content. With initial electric seed field, the dominant effect is the transfer of the electric field energy 
to the hard particle background. The magnetic field energy does increase, albeit non-monotonically, 
without reaching the initial energy density supplied by the seed field. In the case of nonzero initial 
currents, there is a sharp initial increase in the energy density, which is however predominantly 
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FIG. 7: Total energy densities and the contributions from electric and magnetic fields of prolate-phase Weibel 
instabilities with k = O.Itq"^, Tiso/to = 100, and the CGC motivated Debye mass m|, = 1.285/(TisoTo). The 
upper plot corresponds to the initial condition n''(ro) — 1, the lower one to j^{To) ^ 0. The gain rate (times 
To, dotted line) is increased by a factor 5 to make it better visible. 



electric. The magnetic energy density eventually increases, too, again non-monotonically, and only 
slightly higher than in the case with seed electric field. 

We thus find that for our CGC-inspired mass parameter, the Weibel instabilities in the prolate 
phase (where the amount of prolate anisotropy decreases rapidly) are rather weak when compared 
with the required energy densities in the initial seed configuration. 



C. General wave vectors and electric instabilities 

Up to now we have only considered the special situations where the integro-differential equations 
for and were decoupled. We now turn to the more general case where this is no longer the 
case because both k and v are non-zero. Such Fourier components correspond to a physical wave 
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vector whose angle with the r] (or z) axis increases with time according to "& = arctan(rA;/z^). 

In the oblate anisotropic case the growth rate of unstable modes decreases with increasing ^. 
For the a mode the rate tends to zero as i9 — ?• tt/2 and for the "— " mode already at < 7r/4. It is 
therefore of some interest to future 3+1-dimensional simulations of non-Abelian plasma instabilities 
to determine the range of wave numbers k for which the corresponding modes play an important 
part in the evolution. 

Another interesting aspect of the generic case is that it allows us to study electric instabilities 
which appear in the "— " mode. 



1. Electric instabilities for oblate anisotropy 

When solving the coupled integro-differential equations for and for nonzero k and u we 
have again to take care of the nontrivial Gauss law constraint ()3.24p at tq. To do so, we adopt the 
initial data for W field of Eq. (j4.ip also for nonzero k, 

W^{k, vA,y- m) = C5 tanh(2/ - (4.14) 

where the constant C5 is proportional to the initial electric field component parallel to the wave 
vector, 

kU^ + uJr , , 

E\\ = (4.15) 

The electric field, lying in the 1-r] plane, has in general also a component transverse to the wave 
vector, which we shall denote E± in this subsection. (The magnetic field is of course purely 
transverse and points in the 2-direction.) 

An initial W field given by (j4.14p gives rise to an initial charge density (tq) as required by the 
Gauss law constraint, but zero initial spatial currents. In order to also have nonzero initial current 
components and j^, we additionally add the components 

W^^,{k,u;(l),y-r,o) = Ce,7. (4.16) 

As discussed in Sect. IlllCt we can expect to find an electric instability for wave vectors with 
nonvanishing ■d < 7r/4. Since 'd increases with time, we choose small initial values of i?, which over 
some time satisfy the criterion of being within a 45° cone about the r] axis. 

In Fig. [8] the time evolution of the energy density for a mode with nonzero initial and /ctq = 1, 
u = 10 and is shown, for which we can expect an electric instability only for times smaller than IOtq. 
Again we consider a larger mass parameter = 10/(risoTo) to find more significant results and 
indeed we notice that the total energy density rises initially with an significant electric component 
that is almost entirely longitudinal. After the maximum at about 4ro we observe a strong decay 
and only small plasma oscillations after IOtq. 



2. General magnetic instabilities for oblate anisotropy 

For general wave vector and oblate anisotropy, the purely magnetic instabilities reside in the 
modes for all < vr/2, but with vanishing growth rates as — t- 7r/2. Nonvanishing initial 

currents can be simply taken into account by choosing a constant W2{k,v;(j),y — ??o) = Cs, with 

Cg proportional to j^(To). 

Using the same parameters as in Fig. [HI but now with nonzero initial j^, we find an instability 

that is operative up to the somewhat larger time of about 6.5 tq (see Fig. [9]) and which is almost 
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FIG. 8: Total energy density and the contributions from electric and magnetic fields for a mixed A^'^ mode 
with kro = 1 and ly = 10 exhibiting an electric instability. The electric energy density is separated into 
a (dominant) longitudinal and a (small) transverse part with respect to the wave vector. The remaining 
parameters are m|j = 10/(risoTo) and Tiso/tq — 0.1; the initial condition is j^{To; k, v) ^ 0. 
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FIG. 9: Total energy density and its electric and magnetic contributions for the purely transverse mode 
with parameters as in Fig. [8] and initial condition j^(To; fc, v) ^ 0. 



completely in magnetic fields. For this set of parameters it was in fact crucial to have nonzero initial 
currents — with vanishing initial currents and only initial gauge fields we only found decreasing 
solutions. 

Returning to our standard choice of mass parameter, Fig. [10] displays the time evolution of the 
unstable modes with different values of k, the full lines corresponding to the modes, and the 
dashed ones to mixed A^''^ ones. This shows that the most efficient plasma instabilities in the phase 
of oblate anisotropics are concentrated in the range k < 0.2 Tq"^. 
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FIG. 10: Time evolution of unstable modes with — 10 and various values of k (in units of Tq^) and 
Tiso/to ~ 0.1. Full lines correspond to the purely transverse modes, dashed lines to the mixed A^'^ 
modes, both with initial conditions of nonzero initial currents. 



V. SUMMARY AND CONCLUSIONS 

In this paper we have generahzed the semi-analytical analysis of plasma instabilities in an 



anisotropically expanding plasma of Ref. 39(] to general orientations of wave vectors and all possible 
polarizations of the individual Fourier modes. Moreover we have generalized to arbitrary initial data 
in both the collective gauge fields and the W fields of the hard loop formalism, which correspond 
to colored fluctuations in the hard particle distribution and thus directly to the induced currents. 

Besides the well-studied magnetic (Weibel) instabilities of a plasma with oblate anisotropy, we 
have also considered plasma instabilities involving growing electric fields parallel to the wave vector. 
For the latter, the wave vector needs to have a nonzero angle with respect to the axis of expansion, 
which then increases with time, eventually shutting off such instabilities in the expanding case. We 
have also considered a plasma that starts with prolate anisotropy which after some time turns into 
an oblate one. Such instabilities have occasionally been conjectured to be the most interesting for 
isotropization in heavy- ion collisions |12l] . However we found that this type of Weibel instabilities 
requires plasma densities much larger than those suggested by CGC calculations. 



For Weibel instabilities in the oblate phase, Ref. [39|] has previously observed an uncomfortably 
long delay before they overcome the depletion of the energy in initial fields due to the (free- 
streaming) expansion of the plasma. With parameters taken over from CGC calculations, there 
seemed to be very little room for plasma instabilities for the available energies and plasma lifetimes 
at RHIC, while energies and plasma lifetimes expected for heavy-ion collisions at the LHC would 
make an important role conceivable, if the quark-gluon matter to be produced there turns out to 
be sufficiently weakly coupled to behave as a plasma. 

However, Ref. [39| has considered only initial fluctuations in collective fields as seeds for plasma 
instabilities, while physically one should expect fluctuations in both collective fields and in the 
initial hard particle distribution. With our more general initial conditions that allow also for 
fluctuations in the initial hard particle distribution and the corresponding induced currents, we 
find much more favorable conditions for plasma instabilities. When the initial fluctuations in (only) 
induced currents are such that they give rise to the same energy content in collective modes as 
considered in the case of only initial field fluctuations, a surprisingly stark reduction of the delays 
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of plasma instabilities by almost an order of magnitude was obtained. Because of linearity in the 
weak field regime, this implies that the generic case of fluctuations of comparable strength in both 
induced currents and collective fields is overwhelmingly dominated by the modes corresponding to 
only initial current fluctuations. 

In Sect. IIV A 2l we have concluded that for LHC energies the time scales for plasma instabilities 
to set in are of the order of ~ 0.5 fm/c when initial current fluctuations are considered, while for 
RHIC energies these values would be about 2-3 times larger. Although for a significant backreaction 
of the plasma instabilities on the anisotropic hard particle distribution one would presumably have 
to consider times that are somewhat larger, this still seems to keep this mechanism very interesting 
at least for LHC energies. 

Finally, we should emphasize that the present analysis was carried out in the weak-field (linear- 
response) regime. The nonlinear regime of quark-gluon plasma instabilities in the case of boost- 
invariant expansion was considered for an effectively 1+1-dimensional evolution in Ref. [13]. These 
results remain valid as far as the specific non-Abelian dynamics is concerned, but the uncom- 
fortably long delay of the onset of the instabilities largely disappears by considering also initial 
fluctuations in the induced currents (equivalently, in the W fields). However, ultrarelativistic 
plasma instabilities have been found to behave very differently in the regime of nonperturbatively 
large non-Abelian gauge field when a full 3+1-dimensional situation is considered. Work in this 
direction is in progress, for which the present semi-analytical results will provide important cross- 
checks for the initial stages of the evolution of non-Abelian plasma instabilities. 
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Appendix A: Analytical late-time behavior 

For modes with wave vector parallel to the anisotropy direction, i.e. k = and v ^ 0, which 
is the case studied before in Ref. |39j], the expressions for the induced currents simplify, and it is 
possible to study the late-time behavior of single modes analytically by expanding the contributions 
to the memory integrals around r' = r. Late-time behavior in our free-streaming approximation 
means extreme anisotropy, characterized by t^so/t = 0-^1. In this appendix we recapitulate the 
analytical results of Ref. |39!], filling in some details and also show that the late-time behavior is 
not modified by the necessity of including initial values for the longitudinal current in the case of 
longitudinal modes. 

1. Transverse modes 

With A; = 0, the induced currents j^{T]u) are given by Eq. (j4.9p . where for simplicity we set 
-^^{jo'ii^) = and jg(r;i^) = 0. Expanding the integrand of the memory integral around t' = t 
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yields 
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for = Tiso/T ^ 1, the transverse current reduces to 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 



where n = ml^vrriso/S and only terms up to linear order in Tiso/T have been kept. The equation of 
motion for the transverse gauge fields are 



(^drrdr+'^)A\T-u)=f{T-u) 



(A6) 



and acting with d'^r'^ on it we eventually obtain an ordinary differential equation for each mode v 

2 

(d^TdrTdr + u'^d^ + fid^T - ^) A\t; v) ~ 0. 



(A7) 

Simple results for the gauge fields are only obtained for very infrared modes <C 1, where all terms 
proportional to can be neglected, or for high momentum modes ^ 1, where only those terms 
proportional to contribute. We find 



A\t; !/<!)- ci Jo(2V^) + C2yo(2^^S?), 



(A8) 
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FIG. 11: The contribution from to the current proportional to the initial conjugate momentum Jq is negligible 
compared to the rest at late times. This data is for i' = 10 and Tiso/tq = 0.01. 



which is a stable oscillatory solution (Jn{x) and Yn{x) are Bessel functions of the first and second 
kind, respectively), and 

i'(r; 1/ > 1) ~ ci^h{2^) + C2^/fi^i(2V7^), (A9) 

with ci^2 being constants. The modified Bessel functions Kn and /„ have the asymptotic behavior 
Kn{x) ~ exp(— x)/\/2vrj; and In{x) ^ exp(a;)/\/27rx, where the latter describes a rapidly growing 
mode. Therefore we expect that large u modes will be dominant at sufficiently late times with a 
behavior of 

A\t) ~ r^/^ exp(2V7^). (AlO) 

For 1/ ~ 1 the solutions can be written in terms of generalized hypergeometric functions 2-P3 and a 
Meijer G function 39|]. The dominant contribution is 

'3-\/l + 4z/2 3 + Vl + 4z/2 



2 ' 2 

which is compared with the full semi-analytical result in Fig. [2j 



;2,2 - iv,2 + iiy,-fiT), 



(All) 



2. Longitudinal modes 

For longitudinally polarized gauge fields we proceed analogously. In this case it is not a priori 
admissible to drop the the term Jq (r; i^) term in Eq. (j3.20p because of the Gauss law constraint. 
However, numerically, we notice that at late times this contribution is negligible compared to the 
rest of the current, as can be seen in figure [TTl 

Omitting both jQ^T-jv) and the term proportional to j4^(ro;z^), we can approximate Eq. (14. 7p 

by 
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where we have used again ()Aip . but with the factor l/r^^ replaced by 1/r'^. Using 

dy sinh^ y arctan(\/^^^ — 1) 
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and 



+ 6^-2 sij^]^2 -)2 

for 6* = tiso/t ^ 1, we obtain 
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By acting with S^r^ on this expression we find 

dl{r'nT-u)) ^ 2/.a^(M^) _ ^J^A,{r;u), (A16) 



where we can neglect the second part for very large r. The equation of motion for the longitudinal 
gauge fields therefore becomes 

dr^dr + ^)A,iT;u)^0 (A17) 
and the late-time behavior is given by 



T 

This corresponds to stable and oscillatory solutions. 



ci J2(2yV^) + C2Y2i2y^). (A18) 
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